Heavy-Light Decay Constants: Conclusions from the Wilson Action 
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We report on the results of a MILC collaboration calculation of /s , fs s , fn , fn e and their ratios. We discuss 
the most important errors in more detail than we have elsewhere. 



As is well known, precise computations of 
heavy-light decay constants such as fs would 
place stringent constraints on the Standard 
Model. Reference Q describes our recently com- 
pleted evaluation of these decay constant using 
the Wilson action. We have good control of all 
sources of error within the quenched approxima- 
tion. By comparing our quenched results with 
those from lattices with two flavors of staggered 
dynamical quarks, we are also able to estimate 
the error due to quenching. However, for reasons 
described below, the error on this error is proba- 
bly rather large. 

Rather than repeat the full exposition of 
Ref. 0] , we concentrate only on a few key points 
here and discuss them in greater detail than was 
possible previously. This paper should therefore 
be read in conjunction with [Q, which also con- 
tains references to related work. Further, to com- 
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plement the earlier discussion, we focus on our 
data for the ratios of decay constants (/b s / /s and 
rather than the decay constants them- 
selves. 

The largest error within the quenched approx- 
imation comes from (a) the extrapolation to the 
continuum, in conjunction with (b) the chiral ex- 
trapolation, and (c) perturbation theory errors. 
(Errors (a), (b), and (c) cannot be disentanged 
with current data.) Figures [l] and || show /g s / Jb 
and /d,//d, respectively, as functions of lattice 
spacing. In each case, two possible extrapolations 
to the continuum are shown: (1) a linear fit in 
a to all the quenched data, and (2) a constant 
fit to all the quenched data with a<0.5 GeV -1 
(6/<? 2 > 6.0). Since we expect Wilson fermions 
to have 0{a) errors, fit (1) is the most natural 
choice. This choice receives support from the con- 
fidence levels of the fits: in both figures, (1) is a 
much better fit than (2). 
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Figure 1. Results for JbJJb as a function 
of lattice spacing, with central choices for other 
elements of the analysis. The fits are to the 
quenched approximation points (diamonds) only. 



For the decay constants themselves, unlike the 
ratios, the confidence levels of the linear and con- 
stant fits described above are both good, with the 
constant fits in fact having higher confidence lev- 
els than the linear fits. (See Fig. 1 in [jjj.) How- 
ever it would be inconsistent to treat the decay 
constants as independent of a (for a<0.5 GeV -1 ), 
yet fit the ratios linearly. We therefore choose, for 
our central values of both decay constants and 
ratios, the results from the linear fit. Clearly, 
however, the difference between the fits must be 
included in the systematic error estimate. 

As described in Ref . [Q , the errors coming from 
the chiral extrapolations, as well as perturbative 
effects beyond one loop, are entangled with the 
continuum extrapolation error. This is because a 
change in the types of chiral fits used, or a change 
in the one- loop scale q* , moves the individual, 
fixed-a points enough to affect significantly the 
difference between linear and constant continuum 
extrapolations. 

Altogether, we consider 4 choices for the chiral 
fits and 3 choices for q* (see We then com- 
pute each quantity 24 times (2 continuum extrap- 
olations x 4 chiral fits x 3 scale choices) , giving a 
central value and 23 alternatives. The spread in 
the alternatives (taken separately in the positive 
and negative directions) determines the combined 
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Figure 2. 

Same as Fig. [l], but for JdJ Id- 



systematic error due to these three effects. 

For the decay constants, the errors due to each 
of the three effects alone (determined by variation 
in the corresponding choices only) are compara- 
ble, and the combined error computation is non- 
trivial. For example, the combined positive error 
thus determined on fg or fg 3 is ~ 25% smaller 
than the sum in quadrature of the three individ- 
ual errors. This is due to the correlations among 
the errors. 

In contrast, for the ratios of decay constants, 
the errors due to the chiral fits or perturbative 
corrections alone are quite small. This is illus- 
trated for JbJ ' Jb m Fig. |[ Changing the chiral 
fit and/or the q* choice makes very little differ- 
ence, as long as the linear continuum extrapola- 
tion is used in all cases. The combined error for 
each of the ratios then turns out to be, to a good 
approximation, simply the average (over all chi- 
ral and q* choices) of the difference between the 
linear and constant continuum extrapolations. 

Other sources of error within the quenched ap- 
proximation are much smaller than the combined 
error just described and appear to be more or 
less independent of it. They are added to the 
combined error in quadrature. See Ref. [Q for 
details. 

We emphasize that the central values quoted 
come from the quenched approximation. Our 
Np = 2 dynamical fcrmion simulations are used 
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Figure 3. 

Quenched /b,//b with various alternative 
choices in the analysis. 



only to make an estimate of the effects of includ- 
ing virtual quark loops. We feel it would be pre- 
mature to try quote "full QCD" results because: 
(1) the virtual quark mass is fixed and not extrap- 
olated to the chiral limit, (2) we do not believe 
the Np — 2 data is good enough at this point to 
attempt an extrapolation to a — 0, and (3) we 
have two light flavors, not three. 

In practice, the quenching error is estimated 
primarily by comparing the Np = 2 results 
at the smallest available lattice spacing (a ps 
0.45 GeV" 1 , = 5.6, am = 0.01, 16 3 x 32 and 
24 3 x 64) with the quenched results interpolated 
{via the linear fit) to the same lattice spacing. 
The 24 3 x 64 Np — 2 set is new this year and 
gives decay constants which are consistent with, 
but somewhat larger than, the decay constants 
on the older 16 3 x 32 (HEMCGC) set. Since the 
Np — 2 decay constants at this lattice spacing are 
always larger than the corresponding quenched 
ones (see Fig. 1 of Ref. (lj ) , we end up quoting a 
very asymmetric quenching error for f B , /s s , fo , 
and f Ds . For f B JfB and /d„//d, the results 
for Np = 2 and quenched lattices are consistent 
within errors (see e.g., Fig. [l]). It is not surprising 
therefore that the quenching error we estimate in 
this case is therefore roughly the size of the sta- 
tistical errors. 

The effect of changing the way the lattice spac- 



ing is determined (from fixing f n to fixing m p ) 
can also be used as a rough estimate of the size of 
(some) quenching errors. For the decay constants, 
this method gives an error much smaller than the 
direct comparison described above. For f B J f B 
and /d s //d, the two approaches give compara- 
ble estimates (compare Fig. ||with Fig. p. Simi- 
lar statements apply to the effect of changing the 
way the strange quark hopping parameter is de- 
termined (from fixing the K mass to fixing the <fi 
mass). We take the final quenching error as the 
largest error found with any of the methods. 

Our results are (in MeV): 
f B = W7(ll)(tf)(tf); / B3 = l7l(l0)(± 34 )(±f) 
/ Z3 = 192(ll)(t 16 )(t 15 ); / D . = 210(9)(±| B )(±i 7 )- 

The errors are statistical, systematic (within 
the quenched approximation), and systematic (of 
quenching) , respectively. It must be kept in mind 
that the error on the quenching error is large. For 
example, an extrapolation of the Np = 2 results 
to the continuum looks like it would significantly 
increase the (positive) quenching error on the de- 
cay constants. Further, even within the current 
algorithm for estimating quenching, the error on 
the quenching error is >50%, due to the statisti- 
cal fluctuations and the various systematic varia- 
tions in the analysis. 
For the ratios, we find: 

f B Jf B = l.ll(2)(±f)(3) 

fnjfv = l-10(2)(±f)(±i) 

fs/fn. = 0.75(3)(±|)(± 8 ) 

f B Jf Ds = 0.85(3)(t 5 3 )(tJ). 

We note that as experimental measurements of 
fo B improve, the ratios /b//d, and f B J f Da may 
ultimately provide the best way to determine f B 
and f Bs . 
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